CURRENT STATUS TOWARD THE PROTON MASS CALCULATION 

IN LATTICE QCD 

Yoshinobu Kuramashi 
Graduate School of Pure and Applied Sciences and Center for Computational Sciences, 
University of Tsukuha, Tsukuha, Ibaraki 305-8571, Japan 

The proton mass calculation is still a tough challenge for lattice QCD. We discuss the current 
status and difficulties based on the recent PACS-CS results for the hadron spectrum in 2+1 
flavor QCD. 



1 Introduction 

The proton mass calculation has a profound meaning in lattice QCD: To distinguish the proton 
mass from the neutron one we need to incorporate the isospin breaking effects with the different 
up and down quark masses and the electromagnetic interactions. This is still out of reach for 
current lattice QCD calculations. The accomplishment of the first principle calculation of the 
proton mass inevitably means that other physical quantities should be calculated with similar 
precision on the same configurations. 

In this report we show the recent progress in lattice QCD based on the 2+1 flavor lattice 
QCD results obtained by the PACS-CS Collaboration who are currently aiming at the physical 
point simulation^ We discuss the difficulties in lattice QCD from a view point of the systematic 
errors. The toughest problem is the rapid increase of computational cost with the up-down (ud) 
quark mass reduced toward the physical value. We explain why the direct simulation on the 
physical point is required in order to avoid the problems in the chiral extrapolation method. 



2 Difficulties in Lattice QCD Calculation 



Most fundamental quantities in lattice QCD are Green functions in the path-integral formalism: 
(0[C/,g,g]) = ^lvUVqVqO[U,q,q]e-'^culU,<iM^ (1) 

where Sqqy) represents the QCD action defined on the discretized four-dimensional space time. 
U is the so-called link variable which contains the gauge fields, q and q denote the quark and the 
anti-quark fields. Only the Monte Carlo method makes feasible the nonperturbative evaluation 
of the above expression. With appropriate choices of the operator O we can extract various 
physical quantities, e.g., the hadron spectrum. 

There exists two types of errors in lattice QCD: One is the statistical one due to the Monte 
Carlo technique. The other is the systematic ones. The former is arbitrarily reduced according 
to l/\/iV with N the number of the independent configurations (Monte Carlo samples). The 



troublesome is the latter. We have four major systematic errors: (i) finite volume effects, (ii) 
finite lattice spacing effects, (iii) quench approximation and (iv) chiral extrapolation. It is rather 
straightforward to diminish the first and the second errors with the use of larger and finer lattices. 
For almost twenty years after the first lattice QCD calculation of the hadron masses in 1981^ 
most of the large-scale simulations were carried out in the quenched approximation where the sea 
quark effects are neglected. The primary reason is that the 2+1 flavor lattice QCD simulation 
requires O(IO^) times as much computational cost as the quenched approximation. In late 90s 
the CP-PACS collaboration performed a detailed investigation of the quenching effectsP The 
systematic study of the hadron spectrum in the quenched approximation with other systematic 
errors under control reveals that the results deviate from the experimental values at a 10% level. 
The comparison are depicted in Fig. [H where the physical inputs are a set of m^r, mp, niK (closed 
triangles) or mTr,mp,m^ (open triangles) to determine the averaged up-down quark mass, the 
strange one and the lattice spacing a. The confirmation of the discrepancy between the quenched 
results and the experimental values drove us to embark on the 2+1 flavor QCD simulations. 
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Figure 1: Quenched light hadron spectrum com- 
pared with experiment. 



Figure 2: Simulation cost as a function of rriTy/mp. 
See text for details. 



Now the remaining task is to remove the systematic error associated with the chiral extrap- 
olation. Figure [5] illustrates the difficulty: The solid line represents the empirical cost estimate 
for the 2+1 flavor lattice QCD simulation with the Hybrid Monte Carlo (HMC) algorithm 
given by Ukawa in 200131 The cost seems to almost diverge as the m-^jmp ratio approaches the 
physical point. It was obvious that we definitely need not only the increase of the computa- 
tional power but also the algorithmic improvements. Years later the difficulty is overcome by 
the Domain-Decomposed Hybrid Monte Carlo (DDHMC) algorithm.^ Blue circles denote the 
measured computational cost in our simulation armored with several other algorithmic improve- 
ments, which clearly shows that the direct simulation at the physical point is allowed with the 
current computational resources. 

Before physical point simulations, we should examine the logarithmic quark mass dependence 
in the pseudoscalar meson sector predicted by the Chiral Perturbation Theory (ChPT). This 
is a good testing ground to check whether or not the light quark simulations are properly 
performed. In Fig. [3] we plot the ratio w^jra^^ as a function of in lattice unit together with 
the previous CP-PACS/ JLQCD results for comparison. The curvature observed near the chiral 
limit is explained by the SU(2) ChPT prediction: 
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(2) 



where B, f, are the low energy constants and fi is the renormalization scale. Figured compares 



our results for /3j_which is defined by I3 = —64-71^/3 at fi = rriT^, with currentlvavailable data given 
by other groups^A^ Black symbol denotes the phenomenological estimateP Red closed (open) 
symbols are for the results obtained by the SU(3) (SU(2)) ChPT fit on 2+1 flavor dynamical 
configurations. All the results for ^3 reside between 3.0 and 3.5, except for the MILC result 
which is sizably smaller and marginally consistent with others within a large error. 
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Figure 3: m^/niud as a function of mmj. Solid line 
is just for guiding your eyes. 
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Figure 4: Comparison of h obtained by the 2+1 
flavor dynamical simulations. See text for details. 



3 Why is the Physical Point Simulation Necessary? 

Chiral extrapolation with the use of ChPT as a guiding principle is current most popular strat- 
egy to estimate the results at the physical point. The simulation points are usually ranging from 
200 — 300 MeV to 600 — 700 MeV for ttZtt. There are several problems in this procedure. Firstly, 
it is numerically difficult to trace the logarithmic quark mass dependence of the physical quanti- 
ties predicted by ChPT. High precision measurements are required for the reliable extrapolation. 
Secondly, it is not always possible to resort to the ChPT analyses. A typical example is SU(3) 
Heavy Baryon ChPT which completely fails to describe the lattice results for the octet baryon 
massesP Figure [5] shows the next-to-leading order (NLO) fit result for the nucleon mass. This 
difficulty may be practically avoided by the use of the polynomial fit function instead of ChPT. 
We apply a simple linear function of = itiq + arriud + Pi^^s to the lattice data obtained at 
156 MeV < m-TT < 410 MeV. In Fig. [6] we compare our results for the hadron spectrum with 
the experimental values. Most of them are consistent within the error bars, though some cases 
show 2 — 3% deviations at most. Note that we are left with the 0((Aqcd ■ o)^) finite lattice 
spacing effects thanks to the nonperturbative 0(a)-improvement employed in our formulation^ 
This encouraging result, however, does not mean the polynomial extrapolation is a sufficient 
solution. Since we know that mud = is a singular point in ChPT, the convergence radius of 
the analytic expansion around the physical ud quark mass is just < mud < 2m^|j^^"^^', which 
roughly corresponds to MeV < m.^ < 190 MeV. Thirdly, it is impossible to make a proper 
treatment of resonances, e.g. p meson, with the extrapolation method. The reason is quite 
simple: Lattice QCD calculation shows that the pion mass quickly becomes heavier as the ud 
quark mass increases so that the kinematical condition 2m-„: < mp is not satisfied anymore at 
the unphysically large ud quark mass. It is theoretically difficult to predict the real world, where 
the p — > TTvr decay is allowed, by the chiral extrapolation from the virtual world with the decay 

"Similar hadron spectrum is obtained by the BMW Collaboration?^ It is likely that their continuum extrap- 
olation using the simulations at three lattice spacings succeeds in removing the 0(Aqcd • a) errors which are the 
leading finite lattice spacing effects in their formulation. 



forbidden. Fourthly, our final destination is to simulate the different up and down quark masses, 
which is an essential ingredient for the proton mass calculation. The isospin breaking effects are 
so tiny that the reliable evaluation would be difficult with the chiral extrapolation method. 
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Figure 5: 
ChPT fit 
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Figure 6: Liglit liadron spectrum in 2+1 flavor 
QCD with m.^, niK, ma as physical inputs. 



4 Conclusion and Future Plan 

In the history of lattice QCD simulation the computational cost has been a big issue preventing 
us from the direct simulation at the physical point. The chiral extrapolation method, which is 
just a compromise due to the lack of computational power, has intrinsic problems to be avoided. 
Thanks to the rapid increase of computational power and the algorithmic improvements in last 
decade the PACS-CS Collaboration are now able to simulate the physical point directly on a (6 
fm)^. We are going to incorporate the up-down quark mass difference and the electromagnetic 
interactions in the next step. 

Acknowledgments 

Numerical calculations for the present work have been carried out on the PACS-CS computer 
under the "Interdisciplinary Computational Science Program" of Center for Computational Sci- 
ences, University of Tsukuba. This work is supported in part by Grants-in-Aid for Scientific Re- 
search from the Ministry of Education, Culture, Sports, Science and Technology (Nos. 16740147, 
17340066, 18104005, 18540250, 18740130, 19740134, 20340047, 20540248, 20740123, 20740139 ). 

References 

1. PACS-CS Collaboration, S. Aoki et al, Phys. Rev. D 79, 034503 (2009) 

2. H. Hamber and G. Parisi, Phys. Rev. Lett. 47, 1792 (1981). 

3. CP-PACS Collaboration, S. Aoki et al, Phys. Rev. Lett. 84, 238 (2000); Phys. Rev. D 
67, 034503 (2003). 

4. A. Ukawa, Nud. Phys. B (Proc. Suppl.) 106, 195 (2002). 

5. M. Liischer, J. High Energy Phys. 05, 052 (2003); Comput. Phys. Commun. 165, 199 
(2005). 

6. RBC and UKQCD Collaborations, C. Allton et al, Phys. Rev. D 78, 114509 (2008). 

7. C. Bernard et al, Proc. Sci. LAT2007, 090 (2007). 

8. J. Gasser and H. Leutwyler, Ann. of Phys. 158, 142 (1984). 

9. PACS-CS Collaboration, K.-I. Ishikawa et a/., |ar2£iYi0905.0962 [hep-lat]. 



10. Z. Fodor, these proceedings. 



